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Abstract — Centralized automatic train control algorithm is proposed which allows effective adjustment of each train ’s 
departure and running times or intervals when there are unscheduled delays. The novelty lies in that it employs three sets of 
methods, called regulative characteristics, which contain relationships between regulation restrictions ( minimal station dwell 
times and train running times). These relationships are obtained a priori from the simulation of train circulation on the 
detailed railway or metro line model, saving later onboard real-time processing requirements, which is especially useful for 
those train control systems which are being only partially modernized. In the traditional algorithms these regulation 
restrictions are frequently considered fixed, while the changes of system ’s state do affect them, so the usage of the values 
recurrently updated by the regulative characteristics allows to have more flexible running and station times for every train on 
the line. The efficiency of the new algorithm is additionally increased by the usage of forecasting of station dwell times and 
train running times based on the history of these values for each train or station. Proposed algorithm decreases the quantity of 
undesired stops between stations, allowing to save energetic resources spent on repeated train acceleration and to raise, in 
general, the quality of the transportation process algorithm has been tested on Mexico Metro Line simulator confirming 
efficiency increase of 3% for the worst case scenario and up to 10% for the best cases (busiest state of the line plus additional 
temporal speed limits). 

Keywords — Railway transportation, Automatic Train Control, Metro, Regulation Algorithms, Centralized Control, 
Forecasting . 


I. Introduction 

Centralized Automated Train Control (ATC) together with Automatic Train Operation (ATO) systems allow modern 
railways to increment traffic density in order to satisfy the raising passengers affluence while maintaining the desired quality 
of the transportation process in urban and suburban rapid railways. While the problem of rising passenger affluence in all 
kind of urban and suburban transport is being addressed from different points of view [1-7, 12-22] over the years, many 
efforts were spent on the subject of automation of the train control [1, 8-17]. In the core of any ATC-ATO system lies a 
family of algorithms which operates mainly with the values of train running and station well times. Normally [1, 13-25], there 
is at two level hierarchy present: the higher level calculates the timing difference between the desired and executed scenario 
(which is normally based on a schedule, an interval or both). The timing difference is caused by per turbations. If the per 
turbation is considered compensable, the higher level calculates the operating values (running and station times) for each train 
on the line. If the perturbation is non compensable, new scenario (schedule or interval or both) is generated. The lower level 
then executes the control commands generated by the higher level. The way in which the higher level calculates the 
operating times depends on a strategy concerning preferential usage of available timing resources (station time resource 
preference, running time resource preference or hybrid) and on algorithm type: there are algorithms based on schedule which 
operate with astronomical times, and algorithms based on intervals that use relative times, while hybrid algorithms use both 
[1, 10]. A11ATC-ATO algorithms have basic regulation restrictions, such as minimal station well time (or minimal station 
time), minimal interval departure interval and minimal running time. Their absolute minimal values are dictated by line 
and rolling stock physical parameters and traditionally considered as fixed. However, these values can be considered as 
restrictions which depend on current system state, meaning that the changes of departure intervals and running times of all the 
trains in the line caused by perturbations affect them (moving them towards more restrictive values). In this article we present 
the algorithm which considers the modifications in the regulation restrictions caused by system state changes based on 
method, also presented in the article, of regulative characteristics of the station-to- station blocks, where the regulative 
characteristics are families of curves that relate train control parameters and are obtained a priority on the railway line 
simulator. Such approach allows saving processing time for further onboard real-time calculations. Analytical calculations for 
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complete train and line model may involve up to several dozens of differential equations. This kind of calculations is time 
consuming, especially due to the non-linear nature of the line and train parameters. In many cases, experimental modeling is 
recommended over analytical calculation of the functions, especially for many ATC-ATO systems which are being only 
partially upgraded (such as, for example, Mexico City Metro system). We also show the effectiveness of the proposed 
algorithm by simulating its application on a detailed computer model of Mexico City Metro Line and present results with 
several quality criteria of transportation process improvement [1, 23, 25]. Parameters of the line, such as curves, grades, train 
type, blocks lengths and quantity between stations, among others, are considered within the detailed model of the developed 
simulator. 

Proposed algorithm lowers the quantity of undesired stops between stations, allowing saving energetic resources spent on 
repeated train acceleration, for which the saving energy via train regulation is a known strategy [19-28]. The new algorithm 
performs best in schedules with high density and high repeatability, so it might be best suitable for metro lines or relatively 
short suburb lines. However, some of the principles described here might be applicable to longer railway lines. 

II. Traditional Algorithms Examples 

In order to show the importance of considering regulation restrictions as dependable on system state changes, we shall 
analyze two ATC-ATO algorithms: one based on schedule and another based on intervals. 

2.1 Algorithm based on schedule 

One of the simplest algorithms based on schedule compensates train’s arriving time deviation by modifying train’s station 
dwell time, and compensates train’s departure time deviation by modifying station-to-station running time. For example, in 
case of late arrival of the focused [n] train, its real arrival time is greater than the planned: 

& T Aj M = £*/ N - t^[n] >0 


where t^[n] is the real actual time of n train arrival on the j- station; t P -\n\ is the planned arrival time for n-train onto j- 
station; AT a - [ n ] is the deviation of the arrival time for n-train on the j -station. 

Now, the only two timing resources the algorithm has for the train to compensate the late arrival and get back on schedule is 
either to shorten the station dwell time or to shorten the running time to the next station. At first, the algorithm will attempt to 
get back into schedule by shortening only the station time. The train will run at the commanded departure time equal to the 

time planned by schedule t^\ri] — t P D -\n\ in compliance to the regulation restrictions on minimal station time (7^ mm ) and 
minimal departure interval (T^ m ) between the focused train and the ahead train [n-1]: 


tJ/W - l Ajl n - 1] ^ T” 

t C n;-t p Dj [n- l]>T™ n 


( 2 ) 


Dj ujl j ui ( 3 ) 

These restrictions are defined by the safety system of the line, also called Automatic Train Protection (ATP). If both 
restrictions are satisfied, then the timing resource of the station dwell time is enough for getting back on schedule, so, only 
train’s station time is shortened; running time to the next station would not be changed, the perturbation (late arrival) for the 
focused train would be compensated. 


If the condition (2)is violated (calculations demand station time to be shorter than the minimal restriction), it means that the 
timing resource of the station dwell time is not enough to get back on schedule. The generated departure time would be 
ignored as it appears to happen earlier than the minimal station dwell restriction time, thus, the train must depart upon 
fulfilling this minimal station time: 


t C Dj [n] = t R Aj [n] + T r 

t%j[n] = tgj[n - 1 ] + T$ in 


(4) 

(5) 
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At this moment, if either the condition (2) or (3) is violated, then the timing resource of shortening station time is not enough 
for getting back on schedule and the running time resource must be used. So, the station dwell time deviation would be 
compensated by modifying the running timeAT^-[n] on the next station-to-station block: 

A7 % [n] = t p Dj [n] - tg . [n] 

( 6 ) 

while also fulfilling the restriction on minimal running timeT^- : 


r£[n] + Ar£[n]>7- R 7 n (7) 

which is also a minimal restriction defined by the safety system of the line. If the calculated running time is greater than the 
minimal running time, then by executing it the train gets back on schedule by the time of arrival to the next station. 
Otherwise, if the commanded running time appears to be smaller than the minimal, the minimal running time must be 
executed. 

If the timing resource is not enough to return to the planned schedule, then the perturbation is considered non compensable 
and the effective way of controlling the line is by using the algorithm based on intervals (which uses relative times instead of 
astronomical time) [1, 26]. 

2.2 Algorithm based on intervals 

The interval algorithm is effective when the perturbation is greater than the timing resource that the algorithm has for the 
train to get back on schedule. The timing resources that the interval algorithm has for the train delay are the same: shortening 
the station time or shortening the running time. So, if the strict schedule is unreachable than sliding the schedule table ahead 
in time is the way of controlling the line, and the intervals, instead of astronomical timing, can maintain the slid schedule 
together. The general interval algorithm modifies the station time in order to compensate the deviation of arriving 

interval A 7^ [n ] , while it modifies the running time in order to compensate the deviation of departure interval AT^ [n ] : 


AT^M = ^{ATjujlnl&Tfjln- l],A7\f.[n - 2] A T^n-m}} 

ATgj [n] = L 2 {AT DI j[n],AT^j[n- l],Ar a c y [7i- 2] A7^[ji-i]} 


whereL;, L 2 are special operators which depend on the strategy, and m, l - the quantity of a head trains considered in the 
control command generation. For example, in the simplest interval algorithm which considers only one ahead train, the 
departure command is calculated as: 

*Dj M = [«-i]+ T sj M + A 7 s) M (10) 

Where [tt] is the commanded departure time for the n-train, t J - [ft — l] - real departure time of the ahead [n-1] -train, 
r/,= |[n] - planned station time of the n-train, A [ri] - commanded compensation of the station time. 

So the commanded departure time is made up from the sum of real departure time of the last departed train [n-1] added by 
the planned station time of the focused train [n] plus the commanded compensation of the station time modification. All the 
mentioned, once again, while the minimal station time and minimal departure interval are respected (conditions (2) and (3)). 
Otherwise the station and running times to be commanded are defined by the most restrictive of two conditions ((4) and (5)), 
as in the section 2.1. Then, the running time is modified according to (9) for 1 = 1 : 

AT^j = -A T Dlj [n] + k.AT^n - 1], 0 < k, < 1 (l]) 

whereAT^j [n]- is the modification to the running time that is to be commended for the n-train towards the j -station; 
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AT D |j [n] - is the deviation in the departure interval of the n-train (that is waiting to be departed); 

AT^n-l] - is the modification of the running time that was commanded to the ahead train [n-1], and ki is a weight 

coefficient. The calculated modification must alter the commanded running time for the focused train, as long as the minimal 
running time condition (7) is respected. Otherwise minimal running time must be executed. 

As we can see, both examples use fixed regulative restrictions of minimal station time and minimal running time marked by 
the absolute minimum defined by the safety system of the line. However, in the dense traffic of contemporary urban and 
suburban railway systems, and especially metro lines, even little timing perturbations on the line become a cause of raising 
the threshold of these minimal restrictions. For example, the late departure of ahead train would cause a restrictive signaling 
(“red light”) on the arriving of the next running train, if its minimal departure interval and/or its minimal running time values 
are not raised above those absolute minimal values defined by the safety system. In the next paragraph we present a method 
to calculate and use special regulative characteristics of station-to-station blocks, which allows the effective raising of this 
minimum and lowering the probability of restrictive signaling (“red lights”) of the safety system of the line, so the circulation 
without “red lights” appearing for the ahead trains would represent a nonrestrictive circulation. 

III. Introduction of Regulative Characteristics 

In order to consider the way in which the system state changes affect the regulation restrictions, three functions, called 
special regulating characteristics of station-to-station blocks, are proposed. The characteristics are non-linear functions that 
establish the relations (curves) between involved variables of train circulation and are to be obtained a priority during 
experimental simulation process on the line’s model to be later used in the ATO-CTC algorithm database. The values of the 
functions are defined by non-linear aspects of the line: station-to-station block curvature and profile, train acceleration 
characteristics, etc. [1, 12, 15, 17]. Due to this nature, experimental modeling is preferred over analytical calculation of the 
functions. 

The use of the first regulating characteristic will allow nonrestrictive train circulation (no “red lights”) upon the perturbation 
of late departure of the ahead train from the ahead station. In other words, it will prevent a train from being stopped by the 
“red light” between stations if the train ahead is late departing from the station ahead. This situation is typical when the 
passengers behavior on the station delays train’s departure. The regulative characteristic defines the minimal departure 
interval for the focused train [n] as a function of running time and station time of the ahead train [n-1] on the station-to- 
station block ahead: 

= f{T Rj [n - 1 ], T Sj+1 [n - 1 ]} ^ 

where T Rj [n — 1 ] - running time of the ahead train on the station-to-station block ahead, 

T S j+\ \ n ~ 1 ] - stationdwell time of the ahead train on the ahead station. 

The method to calculate one by one the points of the characteristic curve comprehends the following sequence around the 
station-to-station block which is being studied in the simulator: 

1) Process of running and dwelling experimental trains on the simulated block initializes with the absolute minimal 
values of the function defined by the safety system (minimal station dwell time, minimal running time, etc), these 
minimums will be raising their values to form the non-linear function of the regulative characteristic; the station-to- 
station block is limited by the focused station j and ahead station j+1; a virtual depot is created right before station j; a 
sequence of trains are to be commanded out of the depot towards the station j and then towards the station j+1, the 
dwelling and running times of the trains are to be changed during the experiment in order to find the relations that will 
allow nonrestrictive circulation; another virtual depot is created after the ahead station j+1, trains are to be virtually 
disappeared into the second depot, because only one station-to-station block is analyzed at a time, so the circulation 
ahead of j+1 station is not to be studied at this iteration of the process; 

2) First train of the sequence arrives from a virtual depot to the first station of the block, then the train executes the 
minimal station time and departures from the first station executing minimal running time, after arriving to the ahead 
station it executes minimal station time on it and departures towards a second virtual depot, where it disappears; 
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3) While the first train performs its running between station j and station j+1, the second train arrives from the virtual 
depot to the first station and executes minimal departure interval before it departures towards the ahead station with 
minimal (the fastest) running time commanded; now, the second train has two possible outcomes: either it arrives 
successfully to the j+1 station, or it is stopped by a “red light” on the block; 

4) If the second train meets the restrictive signaling(“red light”) on the station-to-station block, then the set of 
commanded values is not suited for nonrestrictive circulation and should not take part of regulative characteristic, it is 
rejected and the process is ready for the next iteration; 

5) The experiment repeats itself with increased value of minimal departure interval until the second and following trains 
(up to a reasonable number, for example, 5 following trains) arrive to the second station without any restrictive 
signaling on the way. At this time one point of one curve of the regulative characteristic is found and saved, i.e. the 
minimal departure interval for the following trains when the ahead trains perform with the current running time and 
the current dwelling time on the ahead station; 

6) At this point the experiment increases the minimal running time by a user-definable step (for example, 5 seconds) and 
the whole sequence repeats itself until the next raised minimal departure interval is found and saved, becoming the 
second point of the first curve of the characteristic; the experiment follows until the whole range of running times is 
studied and the minimal departures intervals for the whole set of running times is found forming a first complete curve 
of the characteristic that corresponds to only one ahead station dwell time; 

7) Now the experiment repeats itself with increased value of ahead station dwell time executing again through the whole 
set of running times, thus, obtaining the second complete curve of the characteristic, point by point; here, the increase 
in the departure time of the ahead train is simulating the perturbation that is causing the train’s delay and the 
maximum value for this delay equals to the maximum perturbation delay that is being studied in the simulation 
experiment; for example, if a dwell time of 5 minutes is selected, then the perturbation delay of 5 minutes is simulated 
and the experiment will run the whole sequence within this range for the matter of ahead station dwell time; 

8) Finally, the complete family of curves is obtained each of them corresponding to a particular ahead station dwell time, 
while relating a range of running times and departure intervals and, thus, a complete first regulative characteristic is 
obtained for the particular station-to-station block which is being studied. Later, the experiment must be repeated for 
each station-to-station block of the line. 

9) The points of the characteristic are obtained for one station-to-station block at a time, but when complete family of 
regulation characteristics obtained for the line it provides the necessary data for the algorithm’s variety which 
considers any number of ahead trains in its calculations. 

“Fig. 1” illustrates the first regulative characteristic graphically (a simplified example with reduced range) for one station-to- 
station block. The family of curves forming the complete first regulative characteristic should be calculated within the widest 
practical range of each involved value that corresponds to the practical perturbation timing. This way the later usage of the 
characteristic in the controlling algorithm will allow circulation without restrictive signaling upon typical perturbations. 


T m [n+1], s 



Fig. I.Example of the first regulative characteristics family for a limited range of station 

TIMES, RUNNING TIMES AND DEPARTURE INTERVALS. 
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While the first regulative characteristic allows nonrestrictive circulation despite departure delays on the ahead station, there is 
other kind of perturbations that should be considered on the line, that is - additional speed limit. Sometimes the additional 
speed limit is introduced when partial segments of the line are being under maintenance. Another reason is during 
precipitations (rain or snow) on the open segments of the line that specially applies to metro lines which are partially located 
underground and partially on the ground, as is the case of Mexico City Metro which was used for main simulation example. 

The use of the second regulative characteristic of station-to-station block will allow non-restrictive circulation upon 
additional speed limit introduction between stations. The characteristic defines the minimal departure interval as a function of 
additional speed limit introduced on the block and station time of the ahead train on the ahead station: 


TSirin] = f{V ADi ,T SJ+1 [n ~ 1 ]} 


(13) 


where V ADj - additional speed limit over station-to-stationblock j, 


Ts j+ An-1] - station time of the ahead train on theahead station. 

The method of calculation follows the similar sequence as with the first regulative characteristic: the minimal departure 
interval is tested for a range of increasing additional speed limits while maintaining constant station time for obtaining one 
curve of the characteristic, then the experiment repeats itself for a range of ahead station dwell times, same range being used 
as in the first regulative characteristic. Later in the algorithm, the minimal departure interval will be selected from the second 
regulative characteristic if there have been an introduction of additional speed limit, or from the first characteristic if no 
additional speed limit is being assigned to the station-to-station block. “Fig. 2” illustrates the example of the reduced range of 
the function graphically. 

T DI [n+l],s 



Fig. 2. Example of the second regulative characteristics family for a limited range of speed 

LIMITS, STATION TIME AND DEPARTURE INTERVALS. 

While the first and second regulative characteristics allow us to effectively select the new minimal departure interval, the 
third characteristic helps to choose the new minimal running time. It defines the minimal running time as a function of 

departure interval T DI [ft] (previously selected either from the first or second characteristic), running time of the ahead 
train T R . [ft — 1] and ahead station dwell time: 


n !?”M = f{T DI jlnlT Bj [n - l],T SJtl [n - 1]} 


(14) 


“Fig. 3” illustrates the function graphically. The method of calculation is similar to the first and second characteristics, but, as 
the third characteristic is a function of three arguments, the simulating sequence requires the whole extra set of iterations for 
the additional argument. The resulting family of curves is three-dimensional, but for the similarity with the first and second 
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characteristics it is presented graphically as several two-dimensional families of curves located side by side on the same 
plane of coordinates (“fig. 3”). 


75 * 85 * 95 

105 115 ■ 125 



Fig. 3. Example of the third regulative characteristics family for a limited range of 

RUNNING TIMES, AHEAD TRAINS RUNNING TIMES, DEPARTURE INTERVALS ORGANIZED IN GROUPS OF 

STATION DWELL TIMES. 

IV. Forecasting station and running times for the ahead trains 

As we can notice from the equations of “section 3”, two arguments of the regulative characteristics require special attention: 
station time of the ahead train on the ahead station T Sj+l [ft — 1] and running time of the ahead train T Rj [ft — 1] . In the 

conditions of dense traffic these values are not always known at the moment of departure of the train for which the regulative 
characteristics are to be selected. Thus, a forecasting of these values is proposed, based on historical values that these 
magnitudes presented. In the proposed model for station time, the deviation of train’s departure is extrapolated upon known 
deviations of real departure times from planned departure times for a series of trains that have already passed the station: 

W n — 1] = ^s/+ — 1] + Fsj+ii n ~ 1] 


where T^ +1 [ft — 1] - planned station dwell time for the ahead train on the ahead station, 

F Sj+ 1 [ft — 1] - forecasted deviation of ahead train’s departure time from the ahead station, 

T S j+\ [ft — 1] - forecasted station time of the ahead train on the ahead station. 

Least square method of first and second order was used for forecasting the deviation based on progression of different series 
of passed trains: 

M 

FJ l+ i[ n - 1] = £ F s R j+1 [n -2-M + iJL 1 

i=0 (16) 

Where F^ +1 [n — 2 — M + i] are the historic values of real deviations for the variable that is being extrapolated; 
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M+l - is the quantity historic trains ahead that are being taken into account for the extrapolation method; 
L 1 - is the operator of least square method of first order (when l = 1) or second order (when 1 = 2). 


There is a necessity of assessing the inaccuracy of the forecasting method, in order to do that we must study the level of 
correlation of the values [l] for a number of historic ahead trains. A set of typical correlations was used: 


Ri(k) = ° 2 


sino) 0 k 

o) 0 k 


R 2 (k) = o 2 e ~ a2k2 
R 3 (k) = cr 2 e -a l k l 


(17) 

(18) 
(19) 


First order extrapolation 




Second order extrapolation 












Fig. 4. Inaccuracy assessment of different correlations with the 1 st and 2 nd order 
EXTRAPOLATIONS: 1 st row 2 nd rowi? 2 ( fe ); 3 rd row i? 3 (fc). 


Where Q~ - is the dispersion of the deviations of F^- +1 [k], so the argument of the correlation is the list of historic ahead 

trains (k=l , 2, . . M+l) and each of them is treated as a moment of time. Assessment of the inaccuracy was performed for 1 st 
and 2 nd order extrapolations (fig. 4), each row contains the comparison of the forecasting accuracy for correlated and non- 
correlated data applied to three different correlation coefficients (Ri(k) = 0.9, 0.8 and 0.7; R 2 (k) = 0.9, 0.8 and 0.7; R 3 (k) = 
0.9, 0.8 and 0.7) that correspond to each of the used correlation bonds. The investigated results here are the values of 
variance (1 st and 3 rd columns) and standard deviation (2 nd and 4 th columns). It can be noted that a number of historic ahead 
trains can be selected for particular correlation coefficient for which the standard deviation of the inaccuracy is lower than 
one, meaning the forecasting is effective. It can be appreciated how fast the standard deviation grows with the decrement of 
correlation coefficient. This allows to define the threshold value of (M+l) historical ahead trains that are to be taken into 
account during extrapolation for which the standard deviation of the forecasted value of [n — 1] is smaller than one. It 

can be seen from the graphics that the threshold of 1 is not crossed for M+l <5. Thus, the uncertainty assessment of the 
method recommends up to 5 ahead trains for effective forecasting of the deviation F. 
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“Fig. 5” shows the block diagram of the forecasting algorithm. Running times of the ahead trains are subjects for fewer 
perturbations as they do not depend on passenger’s behavior on the station, nevertheless their deviations also should be 
forecasted on the same basis in order to effectively satisfy the consideration of regulation restrictions dependence on the 
changes of the system’s state. 

It is worth noting that the forecasting is not the core of the algorithm, but rather an accessory. Other forecasting techniques 
may and should be tested for further improvement of algorithms efficiency. 



Fig. 5. Block diagram of the forecasting algorithm for the values of station times of the 

AHEAD TRAINS 

V. The Complete Algorithm 

The result of adding regulative restriction dependence on system’s state changes as well as forecasting of ahead station and 
running times to the traditional ATO-CTC algorithm is the new proposed algorithm. “Fig. 6” illustrates the block diagram of 
the simple schedule algorithm, discussed in the “Section II” updated with the proposed modifications. The blocks of the 
algorithm make reference to some equations presented earlier in the article. 



Fig. 6. Block diagram of the ATO-CTC algorithm based on schedule updated with regulative 

RESTRICTIONS DEPENDENCE ON SYSTEM’S STATE CHANGES AND THE FORECASTING OF AHEAD TRAINS 

STATION AND RUNNING TIMES. 
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During our investigation other versions of ATO-CTC algorithms were also updated with regulative restrictions dependence 
on system’s state changes and forecasting of ahead trains and running times. This includes algorithms based on intervals, 
hybrid algorithms based on schedule and intervals, and subversions of these algorithms which use different timing resource 
priorities: station time resource priority, running time resource priority, fixed times, station time increasing only, and running 
time increasing only versions of the algorithms. By now, only those algorithms were tested which do not allow running time 
correction “on the go” (after the train has departed the station), because the simulation model resembled a real metro system 
of Mexico City and Moscow (which also resembles many other outdated control systems in the world). Adding the “on the 
go” running time correction possibility is now left for further improvement of presented algorithm and is anticipating further 
improvement of its efficiency. 

VI. Experimental Modeling Results 

New algorithms proposed here were tested on several lines of Moscow Metro simulator and Mexico City Metro simulator, 
also developed by the authors (screenshot presented at “fig. 7”). The selection of these metro systems is explained by the 
origin of the research (Moscow State University of Railway Engineering) from one side and PhD studies of one of the 
authors, supported by Mexican Agency CONACYT, on the other side. Moreover, there are similarities shared by both 
systems: high traffic density, organization of the circulation on the lines, red lights operation. However, there are also 
differences, which are, of course, considered in the detailed model. Mexico City Metro is a system with small number of 
blocks between stations, normally varying from 5 blocks for the closely located stations to 10 blocks for distant ones. 
Another difference is the operation of triple platform terminals in almost all Mexico’s lines. One more diversity is in the fleet 
with rubber wheels for which additional speed limits are applied during rain on open sections. Algorithms were tested with 
different kinds of perturbation and traffic density. A series of quality criteria were used for the comparison of the algorithms, 
the following of them were selected as priority: quantity of emergency stops between stations and specific electrical energy 
consumption of the trains. 



Fig. 7. Fragment of the simulator interface. 

“Fig. 8” illustrates circulation graphic representation of a fragment of the line where single timing perturbation of 20 seconds 
departure delay occurs in one of the stations. Horizontal axe shows the time. Vertical axe shows the stations proportionally 
distanced one from another in the way that correct circulation timing is illustrated as a 45° line. Emergency stops (“red 
lights”) between stations illustrated as dots. As we can see, the response of a traditional algorithm produces three emergency 
stops, while the new version of the algorithm produces only one stop. In case of the new algorithm the first stop is 
unavoidable, as the moment of the first perturbation occurrence is unpredictable. However, new algorithm senses the 
system’s change of state and adjusts its regulative restrictions (by raising the departure interval and/or the running time of the 
next train), which helps to avoid the next two emergency stops. 
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Fig. 8. Executed circulation graphic of the fragment of Mexico City Metro Line 3 with 

TRADITIONAL (LEFT) AND NEW (RIGHT) ALGORITHM USED. 


“Fig. 9” illustrates line’s behavior with traditional and new algorithms upon multiple timing perturbations on several stations 
during 15 minutes period of time. As seen before, new algorithm detects the system’s change of state and manages to avoid 
the majority of emergency stops. In fact, the quantity of emergency stops equals the quantity of trains which already left the 
station and cannot already correct its departure interval to avoid the emergency stop ahead. This means that the reached 
emergency stops quantity is the minimum for an algorithm with no running time correction “on the go”. 


Emergency stops quantity reduction decrements the times train has to accelerate on the station -to -station block. “Fig. 10” 
shows the specific energy consumption rates upon repetitive timing perturbations for a variety of eight traditional and eight 
new algorithms (marked with *). As we can see, the new family of algorithms shows noticeable consumption reduction. In 
fact, energy savings reached from 3 to 10%, depending on the kind of perturbation, the density of the traffic and the 
algorithm tuning used. 
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Fig. 9. Executed circulation graphic of the fragment of Mexico City Metro Line 3 with 

TRADITIONAL (UP) AND NEW (DOWN) ALGORITHM USED. 
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Fig. 10. Specific energy consumption comparison for a family of traditional ATO-CTC 

ALGORITHMS AND NEW PROPOSED ALGORITHMS (MARKED WITH *) BASED ON SCHEDULE (S), INTERVALS (I) 

AND HYBRID VERSIONS. 

VII. Conclusion 

Results and methods conceived are mostly applicable for metro lines and can be kept for suburban railway lines with the 
stations located close to each other, so there is small deviation in running times and there are small differences in service 
patterns so that the dwell time can be forecast. 

The core block of the new algorithm is based on the regulative characteristics, which are calculated a priori on the line 
simulator, while the forecasting only improves its efficiency in the conditions of highly dense traffic (that is when sequential 
perturbations are most likely correlated between each other). It might be noted that the usage of other forecasting techniques 
[2, 22, 28] could be adopted and studied in order to work within regulative characteristics methodology. In any case, 
following the presented methodology, the great majority of calculations are performed only once every time there is a 
constructional update on some station -to -station block, for which the regulative characteristics must be also updated, because 
the simulation parameters would change. This leaves few calculation cycles to be performed “on the go” and makes the new 
algorithm particularly suitable for real-time regulation especially on outdated control systems with no great calculation- 
power onboard, which is the popular case of partial modernization. On the other hand, energy saving and general 
transportation process quality improvement makes presented algorithms recommendable for practical usage. 

Acknowledgements 

This work was supported in part by Consejo Nacional de Ciencia y Tecnologia (CONACYT, Mexico) under grant 210646 
for Pavel Vorobiev’s PhD studies in Moscow State University of Railway Engineering. 

REFERENCES 

[1] L.A. Baranov and V.G. Sidorenko. Subway traffic modeling, P roc. International Scientific and Practical Conference “ The techniques 
ofrailway traffic control at the beginning of 21th century ”, Warsaw, Poland, 1999, pp. 23-26. 

[2] Yangsheng Jiang, Juanxiu Zhu, Lu Hu, Xi Lin and AfaqKhattak.A G/G(n)/C/C state -dependent simulation model for metro station 
corridor width design. Journal of Advanced Transportation 2015. DOI: 10.1002/atr.l315. 

[3] Ya-Juan Chen, Zhi-Chun Li, and William H. K. Lam. Modeling transit technology selection in a linear transportation corridor. 
Journal of Advanced Transportation 2015; 49(1): 48-72. DOI: 10.1002/atr.l262. 

[4] Xiaoming Xu, Keping Li and Xiang Li. A multi-objective subway timetable optimization approach with minimum passenger time and 
energy consumption. Journal of Advanced Transportation 2015. DOI: 10.1002/atr.l317. 

[5] PrabhatShrivastava, S. L. Dhingra and P. J. Gundaliya. Application of genetic algorithm for scheduling and schedule coordination 
problems. Journal of Advanced Transportation 2010; 36(1): 23-41. DOI: 10.1002/atr.5670360103. 

[6] J. Gonzalez, C. Rodriguez, J. Blanquer, J. M. Mera, E. Castellote and R. Santos. Increase of metro line capacity by optimisation of 
track circuit length and location: In a distance to go system. Journal of Advanced Transportation 2010; 44(2): 53-71. 
DOI: 10.1002/atr.l09 

[7] Bin Yu, Zhong-Zhen Yang, Kang Chen and Bo Yu. Hybrid model for prediction of bus arrival times at next station. Journal of 
Advanced Transportation 2010; 44(3): 193-204. 

[8] Asuka, M ; Kataoka, K ; Komaya, K ; Nishida, S . Automatic Train Operation Using Autonomic Prediction of Train Runs , Electrical 
Engineering in Japan, vol. 175, issue 3, pp.65-73, 2011, DOI: 10.1002/eej.21080. 

[9] Zhu Jin Ling , Li Yue Zong , Peng Ling Wang , Qing Yuan Wang . Research of the automatic train operation algorithm. Proc. Of Control 
Conference CCC, pp. 8161 - 8166, 26-28 July 2013. 

[10] SACEM: The Advanced Train Control System, MATRA Transport, Paris, 1985. 


Page | 95 


International Journal of Engineering Research & Science (IJOER) 


ISSN: [2395-6992] 


[Vol-2, Issue-4 April- 2016] 


[11] R.-M. Goodall and W. Kortum. Mechatronic developments for railway vehicles of the future. Control Eng. Pract ., pp. 887-898, Oct. 
2002. 

[12] Hairong Dong, Bin Ning, BaigenCai, ZhongshengHou. Automatic train control system development and simulation for high-speed 
railways, Circuits and Systems Magazine , IEEE, vol. 10, issue 2, pp. 6-18, May 2010. 

[13] S. Yasunobu, S. Miyamoto, and H. Ihara, Fuzzy control for automatic train operation system .Proceedings of Int. Conf. on 
Transportation Systems , pp. 33-39, 1983. 

[14] T. Tang and L.-Y. Huang. A survey of control algorithm for automatic train operation. China J. Railway Soc., vol. 25, issue 2, pp. 98- 
102, 2003. 

[15] H.-W. Lawson , S. Wallin , B. Bryntse and B. Friman. Twenty years of safe train control in Sweden. Proc. IEEE Int. Conf. Computer 
Based Systems , pp.289-293, 2001. 

[16] C. Binard and M.-V. Liefferinge. ATBL0-A first step towards an European vital computer (EYC) for ATC. Proc. IEEE Electric 
Railways in a United Europe, pp. 116-120, 1995. 

[17] M.-K. Banerjee and N.-E. Hoda. Review of the automatic train control system for Cairo Metro Line 2. Power Eng. J., pp. 217 -228, 
1998. 

[18] R.-D. Pascoe and N.-T. Eichom. What is communication-based train control. IEEE Veh. Technol. Mag., pp. 16 -21, 2009. 

[19] V.-A. Minin , V.-A. Shishliakov , J.-N. Holyoak , D.-A. Johnston and N.-A. Lepsky. Development of the communications-based train 
control system for Moscow metro. Proc. ASME/IEEE Railroad Conf, pp.201-210, 1997. 

[20] S.-H. Han , Y.-S. Byen , J.-H. Baek , T.-K. An , S.-G. Lee and H.-J. Park. An optimal automatic train operation (ATO) control using 
genetic algorithms. Proc. IEEE Region 10 Conf. TENCON, vol. 1, pp. 360-362, 1999. 

[21] L. Jia , X. Zhang and Z. Xie. Automatic train control: an intelligent approach. Proc. IEEE Computer, Communication, Control Power 
Engineering, vol. 4, pp. 338-342, 1993. 

[22] Jih-Wen Sheu, Wei-Song Lin. Energy-saving automatic train regulation using dual heuristic programming. IEEE Transactions on 
Vehicular Technology, vol. 61, issue 4, pp. 1503-1514, May 2012. 

[23] W. Gunselmann. Technologies for increased energy efficiency in railway systems. Eur. Conf. Power Electron. Applications , 2005. 

[24] Institute for Future Studies and Technology Assessment. “Evaluation of energy efficiency technologies for rolling stock and train 
operation of railways (EVENT”). Final report, Berlin, 2003. 

[25] R. Anderson. Maximizing the potential for metros to reduce energy consumption and deliver low-carbon transportation in 
cities. MetroRail Asia, 2009. 

[26] S.-C. Chang and Y.-C. Chung. From timetabling to train regulation — A new train operation model. Inf. Softw. Technol., vol. 47, 
issue 9, pp. 575 -585, 2005. 

[27] Y. Ding and S. I. Chien. Improving transit service quality and headway regularity with real-time control. Transp. Res. Rec., J. Transp. 
Res. Board, vol. 1760, pp. 161 -170, 2001. 

[28] A. Fernandez, A. P. Cucala, B. Vitoriano and F. d. Cuadra. Predictive traffic regulation for metro loop lines based on quadratic 
programming. Proc. Inst. Mech. Eng. F, J. Rail Rapid Transit, vol. 220, issue 2, pp.79 -89, 2006. 


Page | 96 




